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When calls offered to a primary group of trunks find all of them 
busy, provisions are often made for these calls to overflow to other 
groups of trunks. Such traffic overflow systems have been of interest 
for a long time, but recently overflow systems that allow for some 
calls to be queued have been of importance. In this paper we analyze 
a traffic overflow system with queuing, which consists of a primary 
and a secondary group. The system which we consider here differs 
from the two systems we investigated earlier, in that no overflow from 
the primary to the secondary is permitted if there is a waiting space 
available in the primary queue. As with the earlier investigations, we 
adopt an analytical approach which considerably reduces the dimen- 
sions of the problem, and simplifies the calculation of various steady- 
state quantities of interest. Our results include expressions for the 
loss probabilities, the average waiting times in the queues, and the 
average number of demands in service in each group. 

I. INTRODUCTION 

In this paper a particular overflow system with queuing is analyzed. 
The system consists of two groups, a primary and a secondary, with n k 
servers and q k waiting spaces, which receive demands from indepen- 
dent Poisson sources S k with arrival rates X* > 0, k = 1 and 2, 
respectively, as depicted in Fig. 1. The service times of the demands 
are independent, and exponentially distributed with mean service rate 
(i > 0. If all n 2 servers in the secondary are busy when a demand from 
S 2 arrives, the demand is queued if one of the q 2 waiting spaces is 
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available, otherwise it is lost (blocked and cleared from the system). 
Demands waiting in the secondary queue enter service (in some 
prescribed order) as servers in the secondary become free. 

If all n i servers in the primary are busy when a demand from S\ 
arrives, the demand is queued in the primary, if one of the qi waiting 
spaces is available. No overflow is permitted from the primary queue, 
so that a demand in the primary queue must wait for a server in the 
primary to become free. If all rti servers in the primary are busy and 
all q\ waiting spaces are occupied, when a demand from S\ arrives, the 
demand is served in the secondary, if there is a free server and there 
are no demands waiting in the secondary queue, otherwise it is lost. 

The overflow system described above differs from the two systems 
which we investigated earlier, 1,2 in that no overflow from the primary 
to the secondary is permitted if there is a waiting space available in 
the primary queue. This restriction was one invoked by Anderson. 3 In 
the two systems investigated earlier, arriving calls can overflow when 
the primary queue is not full. The system considered in this paper is 
a particular case of the one considered by Rath, 4 which was composed 
of two queues, one of which is allowed to overflow to the other under 
specified conditions involving the queue lengths. He obtained some 
numerical solutions using a Gauss-Seidel iteration technique, but none 
of these correspond to the particular system that we are considering. 
He also developed an approximate procedure for analyzing his system, 
based on the use of the Interrupted Poisson Process. 

Here we analyze the overflow system using a technique analogous 
to the one introduced in the earlier paper. 1 Let py denote the steady- 
state probability that there are i demands in the primary and j 
demands in the secondary, either in service or waiting. These proba- 
bilities satisfy a set of generalized birth-and-death equations, which 
take the form of partial difference equations connecting neighboring 
states. We carry out an analysis that reduces the dimensions of the 
problem, which may be considerable in cases of interest. The basic 
technique is to separate variables in the region away from a certain 
boundary of state space. This leads to an eigenvalue problem for the 
separation constant. The eigenvalues are roots of a polynomial equa- 
tion. The probabilities ptj are then represented in terms of the corre- 
sponding eigenfunctions. The constant coefficients in these represen- 
tations are determined from the boundary conditions and the normal- 
ization condition that the sum of the probabilities is unity. 

Various steady-state quantities are of interest, which may be ex- 
pressed in terms of the probabilities />i>. The quantities include the loss 
(or blocking) probabilities, the average waiting times in the queues, 
and the average number of demands in service in each group. These 
quantities may be expressed directly in terms of the constant coeffi- 
cients which occur in the representations for the probabilities py. Thus 

2 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1 981 



the steady-state quantities of interest may be calculated directly, 
without having to calculate the probabilities p,> Here again the reduc- 
tion in the dimensions of the problem is valuable. 

Only the theoretical results are presented in this paper. Numerical 
results will be presented in a forthcoming paper by Kaufman, Seery, 
and Morrison. 5 Results will be given there for the two overflow systems 
considered previously, based on the earlier analysis, 1 as well as for the 
system considered in this paper. 

Section II discusses the representation of the probabilities p„ in 
terms of the eigenfunctions, and the boundary and normalization 
conditions. Various steady-state quantities of interest are calculated in 
Section III. The appendix gives properties of the eigenfunctions. 

We assume throughout the analysis that qi > 1, since the system 
considered in this paper, and the two systems analyzed earlier, are 
identical if q\ — 0, i.e., if there is no primary queue. However, the 
results of this paper reduce to those obtained earlier 1 if qi = 0. If q 2 is 
large, or even infinite, an alternate analysis, analogous to that pre- 
sented for the other two systems, 2 may be carried out for the present 
system, but we do not pursue that here. 

II. REPRESENTATION AND BOUNDARY CONDITIONS 

We let pij denote the steady-state probability that there are i 
demands in the primary and j demands in the secondary, either in 
service or waiting. These probabilities satisfy a set of generalized birth- 
and-death equations, 6 which may be derived in a straightforward 
manner. We define the traffic intensities 

cti = Ai/ju, a 2 = A2/M. (D 

and let the total number of servers and waiting spaces in each group 
be 

ki — Hi + q\, k 2 = n 2 + q 2 . (2) 

It is convenient to introduce the function 

-I 1 ' *** (3) 



x/ \0, l<0, 
as well as the Kronecker delta 

*-{S i'A (4 > 

Then the birth-and-death equations are 
[ai(l - 8 tki xj-n 2 ) + a 2 (l - 6> 2 ) + min(i, /ii) + min(y, n 2 )]p,j 
= ai(l - 8>o)pi-ij + (1 - SjoHaiSikjfa-j + a 2 )pi,j-\ 
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+ (1 - 8ik,)mm(i + 1, nfipMj + (1 - S jk2 )mm(j + 1, n 2 )pu+h (5) 
for < i < Ai, <y < A 2 . The normalization condition is 

ZSA/-L (6) 

i=0 /=0 

For i ^ Ai, the variables in (5) may be separated, and there are 
solutions of the form aiflj, where 

[d + min(i, Hi) + p]ai = ai(l - S,o)a.-i + min(c + 1, /ii)a,+i, (7) 

for < i < Aii - 1, and 

[a 2 (l - 5/a 2 ) + mint/, n 2 ) - p]ft 

= a 2 (l - «/))A--i + (1 " WminO" + 1, n 2 )/?, +1 , (8) 

for <y < A 2 , and p is a separation constant. Hence, from (7), 

(ai + i + p)ai = a a (l - 5a))a,-i + (i + l)a,+i, (9) 

for < i < rii — 1. The solution of (9) may be expressed in terms of 
Poisson-Charlier 7,8 polynomials. We here denote the solution of (9) for 
which oto = 1 by Si(p, a\). The properties of Si(p, a) which we will need 
are given in the appendix. 
We assume that (ft > 1. Then, from (7), 

(ai + n x + p)oti = aioii-i + nm+i, (10) 

for /ii < i < k\ — 1. The solution of (10) may be expressed in terms of 
Chebyshev polynomials of the second kind, 9 Ui(x). It is convenient to 
define 

The appendix gives the properties of these functions that we need. We 
note here, however, that U {x) = 1 and U-i(x) = 0. From (9), (10), 
(53), and (64), with a suitable normalization, it follows that 



<Xi(p) = 



— I Si{p, ai), < i < Hi, 

— J [s ni (p, ai)S2,- ni (p)-s ni -i(p, ai)Q,-- ni -i(p)], 



(12) 



ni < i < *i. 
Next, from (8), 

(a 2 +> - p)A- = a 2 (l - Mft-i + 1/ + DAVi, (13) 

for < y < n 2 — 1. It follows from (53) that fij is proportional to 
Sj{-p, ai) for £/£ rii. Also, 
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[a 2 (l - 8jk 2 ) + n 2 - p]fo = a-£j-i + n 2 (l - 9/hJPj+u < 14 ) 

for n 2 <,/ < k 2 . Corresponding to (11), we define 

We also define 

hip) = ¥*,-» - *v*-i(p). (16) 

The appendix gives the properties of these functions that we need. It 
follows from (14) and (62) that ft is proportional to <fe(p) for n 2 - 1 < 
7 < k 2 . 

Consequently, we take 



O t x _ f M-P. A2)<>n 2 (p), <> < /12, 

Wl Is^-P, aa)*-(p), n 2 - 1 < y < fe, 



(17) 



where 

8n8-l(-p, a 2 )4>n^p) = fci^-p, fl2)*nr-l(p)- ( 18 ) 

This equation may be written in the form 

p[s„ 2 (l - p, a 2 )% 2 (p) - s^-id - P, 02)**-i(p)] = 0. (19) 

The expression in the square brackets in (19) is a polynomial in p of 
degree k 2 = n 2 + q 2 . It was shown 1 that its zeros are positive and 
distinct, and we denote them by p m , m = 1, • • • , k 2 . We also let p = 0. 
It follows that we may represent the probabilities />,, in the form 

2 C m Cii(pm)Sj(-p m , a^n^pm), <> < rt 2 , 

m-0 

(20) 

*2 

X Cmtt^pmjs^t-Pm, a 2 )(f>j{p m ), Tl 2 < > < A 2 , 



Pv = 



for < i < Ai, where a t (p) is defined in (12), and the constants c m are 
to be determined. 

It remains to satisfy the boundary conditions corresponding to i = 
ki in (5), as well as the normalization condition (6). If we set i = ki in 
(5), we obtain 

(d + a 2 + /ii +J)pk x j = OiP*,-i.> + («i + a z)(l - 8jo)Pk u j-i 

+ (j + l)Pk t .j+i, (21) 

for < 7 < n 2 - 1, 

[a 2 (l - 5g 2 ,o) + Ml + «2]pA,.n, 

= a\Pk { -\,n, + (ai + a 2 )p* 1 ^ 2 -i + « 2 (1 - 5 92 .o)p* 1 ^ 2 +i, (22) 
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and, if q 2 > 1, 

[a 2 (l - 8jk 2 ) + Mi + n 2 ]pk x j 

= aipk t -i,j + a 2 pk i ,j-i + n 2 (l — &jk)pk v j+u (23) 

for n 2 + 1 < 7 < A 2 . 

If we substitute (20) into (21), after reduction with the help of the 
recurrence relations in the appendix, we find that 

X c m {p m [s„,(l + p m , ai)J2 7l (p m ) - s„j-i(l + p m , ai)fl 9l -i(p m )]s/(-p m , a 2 ) 

m-0 

+ ai[s„,(p m , ai)£L Ql {p m ) - s ni -i(p m , ai)S2 9l -i(p m )]s,(-l - p m , a 2 )} 

•**,(/»») = 0, (24) 
for <y < /i 2 - 1. Also, from (23), it is found that 

*2 

^ Cmp m s^,(-p m , a 2 )[s„,(l + pm, ai)fl 9l (p m ) 

m-0 

- s„,-i(l + p m , aiXVitPaHMP*) = 0, (25) 

for n 2 + 1 < y < k 2 . It may be shown that the boundary condition (22) 
is redundant, as is to be expected. Thus the constants c m are deter- 
mined by (24) and (25) only to within a multiplicative constant, which 
is determined from the normalization condition (6). 
From (20), with the help of (16), (19), (57), and (58), it is found that 

*2 

2 PU = co«,(0)[s„ 2 (l, a 2 )* 92 (0) - Snt-id, a 2 )^ 2 -i(0)], (26) 

y-o 



for < i < ki. But, from (12) and (66), 



a,(0) = 



— ) s,(0, d), 



< i < »i, 



tr- 



(27) 



— ) s„,(0, ai), ni < i < Ai. 



Hence, from (26) and (27), with the help of (57), (58), and (65), the 
normalization condition (6) implies that 

co[s„,(l, a,)fl,,(0) - s„,_i(l, ai)SVi(O)] 

• [s^(l, a 2 )^ 2 (0) - s^-id, fb)¥«-i(0)] = 1. (28) 

Once the constants c m have been determined, the steady-state prob- 
abilities pu may be calculated from (20). We remark that the number 
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of constants to be determined is only k 2 + 1, whereas the number of 
probabilities p, y is (k\ + l){k 2 + 1). 

III. SOME STEADY-STATE QUANTITIES 

We proceed now to the calculation of various steady-state quantities 
of interest. These quantities are depicted in the diagram of Fig. 1, 
which indicates the mean flow rates. The loss probabilities L\ and L 2 
are given by 

k, *, 

U = 2 p*,./, L 2 = £ pij* (29) 

j-n. 2 <=0 

and the probabilities that a demand from the primary, or secondary, 
source is queued on arrival are 

Qi = Y i pih q 2 = (i - s 9 ,o) i i pu. oo) 

i-n, y=0 (=0 /-/fcj 

The probability that a demand arriving from the primary source 
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Fig. 1— Mean flow rates for an overflow system with queuing; Poisson arrival rates 
Xi and X2, loss probabilities L\ and L 2 , queuing probabilities Q\ and Q 2 , and overflow 
probability /12. 
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overflows (immediately) is 

/12 = J Pk lt j. (31) 

7=0 

Since the mean service rate is /u, the mean departure rate from the 
primary queue to the primary servers is 

Rn = n lf i £ Ipij, (32) 

i— n,+l >=0 

while the mean departure rate from the secondary queue to the 
secondary servers is 

*1 *2 

R 22 = ntp(l - 6 92 . ) I S Pij. (33) 

1=0 7=^+1 

The average number of demands in the primary and secondary queues 
are 

Vi - J 2 (i - ni)pij, Vi-J JO'- mfc/. (34) 

/— /i j 7—0 i— 7—^2 

Also, the average number of demands in service in the two groups are 

Xx = £ X min(i', nOp.y, X 2 =H mmO'. HaJ/ty. (35) 

1-0 7-0 1=0 7=0 

Now, according to Little's theorem, 6 the average number of demands 
in a queuing system is equal to the average rate of arrival of demands 
to that system times the average time spent in that system. If we apply 
this result to the primary and secondary queues, we find that the 
average waiting times of the demands which are queued in the primary 
or in the secondary are given by 

W 1 =-^-, W 2 = -^- (q 2 >D, (36) 

Ait^i A2*^2 

respectively, independently of the order of service within each queue. 
Also, if we apply Little's theorem to the primary and secondary groups 
of servers, we obtain 

Ml - U - In) = 11X1, A 2 (l - L 2 ) + Ai/i 2 = vX 2 . (37) 

The steady-state quantities of interest may be expressed in terms of 
the constants c m with the help of the representations in (20). From 
(29) it is found, with the help of (12) and (16), that 

*2 

L\= 2 c m [s ni (p m , ai)ft 9l (p m ) - s„,-i(p m , ai)ft 9l -i(p m )] 

m-0 

•s n2 {-p m ,a 2 )% 2 (p m ). (38) 
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We define 

do = c [s„ 2 (l, a 2 )% 2 (0) - Sn 2 -i(l, a 2 )% 2 -i(0)] 

= [s„,(l, 01)0^(0) - Sn.-id, Oi^-iWr 1 , (39) 

from (28). Then, from (30) it is found, with the help of (26), (27), and 
(65), that 

Q x = dos ni (0, Oi)[Q,,(0) - 1]. (40) 

Moreover, from (29) and (31), it follows that 

Li + In = rfos„,(0, oi), (41) 

and from (32) it follows that R u = XiQi, as is to be expected, since in 
the steady state the mean departure rate from the queue is equal to 
the mean arrival rate to it. 
We define 

\q -(q+ 1)| + £* +1 ]/(l - 2 , 1+ 1, 
fcgfo + 1), £=1. K ' 



A,«) = S «*-' = 



9 

9-1 - J 



/=1 

Then, from (34), with the help of (26), (27), and (39), it is found that 
ri-dos Bl (0,ai)A,/^J. (43) 

Also, from (35), with the help of (56), (58), (59), and (65), it follows 
that 

X, = rfoflx (j J s„,-i(l, a,) + s„,(0, a,)[a„(0) - 1] J. (44) 

It may be verified, with the help of (39), (57), and (65), that (41) and 
(44) are consistent with (37). The explicitness of the expressions for 
the quantities in (40), (41), (43), and (44) is due to the fact that these 
quantities are not affected by the secondary. This, of course, is not the 
case for the loss probability L\, which is given by (38). 

Next, from (29), since <f>k 2 (p) ■ 1, it is found, with the help of (20) 
and (68), that 

L 2 = £ c m Sn 2 (-pm, a 2 )[s n ,(l + p m , ai)Sl gi (p m ) 

m=0 

- s n ,-i(l + p m , ai)S2 9l _i(p m )]. (45) 
Also, from (35), with the help of (6) and (59), it follows that 

*2 

X 2 = n 2 ~ £ C m S„ 2 -i(2 - p m , 02)^n 2 (Pm) 

•[S n ,(l + pm, ai)fi 7| (p m ) ~ S„-i(l + p m , Ol)0,,-l(p m )]. (46) 
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In view of (38), (41), (45), and (46), the second relationship in (37) 
provides a useful numerical check. 

We now define 

*i 
rj-^Pij, n 2 <y<A 2 . (47) 

«=0 

If we sum on i in (5), we obtain 

[a 2 (l - 6> 2 ) + /i 2 ]r, = ChTj-i + n 2 (l - 8 jkl )r j+h (48) 

for n 2 + l<j^ k 2 . It follows that 

Haf> = a 2 r y _i, w 2 + 1 <y < k 2 . (49) 

Hence, since L 2 = r* 2 , from (29) and (47), 

r ; -(-] 2 U n 2 ^j<k 2 . (50) 

Then, from (30) and (33), with the help of (63), we obtain 

Q 2 = [% 2 (0) - 1]L 2 , (51) 

and # 22 = \ 2 Q 2 , as is to be expected. Also, from (34) and (42), it follows 
that 

This completes the calculation of expressions for the steady-state 
quantities of interest. 
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APPENDIX 

We define s,(A, a) by the recurrence relation 

(a + i + X)si(\, a) = a(l - 5,o)s,-i(A, a) + (i + l)s« + i(A, a); 

So(K a) = 1, (53) 

for i = 0, 1, • ■ • . Thus s„(A, a) is a polynomial of degree n in both A and 
a, and it may be related to a Poisson-Charlier polynomial. 7,8 However, 
we give here the properties of s„(A, a) which we will need. An explicit 
formula is 1 

^■* = i§^: (54) 
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where 

(X) = 1, (X) r = \(X + 1) • • - (A + r - 1), r = 1, 2, . . • . (55) 
It was also shown 1 that 

(i + l)s I+1 (A, a) = asi(\, a) + Xsi(X + 1, a) (56) 

and 

Si (X, a) = s,(A + 1, a) - (1 - fo)«-ift + 1, a). (57) 

From (57) it follows that 

l Si (X, a) = s n (X + 1, a), (58) 

i=0 

and, from (56) and (58), we deduce that 

£ (n - 3a(A, a) = (1 - 5 n0 )s n -i(A + 2, a). (59) 



1=0 



We now turn our attention to the Chebyshev polynomials of the 
second kind, 9 Ui(x). They may be defined by the recurrence relation 

2xUi(x) = U M (x) + Ui-M; U-Ax) = 0, U (x) ■ 1, (60) 
for Z = 0, 1, • • • . From (15) and (60) it follows that 
(a 2 + n 2 - p)"9i(p) = a 2 ^/+i(p) + n 2 */-i(p), 

^_x(p) = 0, *o(p) - 1. (61) 

From (16) and (61) we deduce that 

[a 2 (l - 8jk 2 ) + n 2 - pWp) = M>,-i(p) + n 2 (l - Sjkjh+iip), (62) 
for j < k 2 . Also, from (61), it may be shown by induction that 

*/<0) = 1 (-Y. (63) 

r=o \a 2 ) 

Next, from (11) and (60) it follows that 

(a\ + n\ + p)tii{p) = a&i+i(p) + nitii-i(p), 

fi-i(p) ■ 0, &o(p) - 1. (64) 

In particular, it is found by induction that 

0/(0) = i(-)- ( 65 > 

Since s,(0, a) = a'/i\, from (54), it follows that 
s„,(0, ai)G/(0) - s„,-i(0, ai)B/_i(0) = s„,(0, a,), / = 0, 1, • ■ • . (66) 
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Next, from (9) and (10), we deduce that 

p £ oti(p) = (m + p)a kl (p) - aia*,-i(p). (67) 

Then, with the help of (12), (56), (57), and (64), it may be shown that 

*! 

I nip) = s„,(l + p, Oi)B,,(p) - s n ,-i(l + p, ai)n„_,(p), (68) 

i-O 

for p^O. Moreover, this result holds for p = also, from continuity. 
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